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Abstract 

The self-consistent theory of the correlation effects in Highly Correlated Sys- 
tems(HCS) is presented. The novel Irreducible Green's Functions(IGF) method is 
discused in detail for the Hubbard model and random Hubbard model. The inter- 
polative solution for the quasiparticle spectrum, which is valid for both the atomic 
and band limit is obtained. The (IGF) method permits to calculate the quasi- 
particle spectra of many-particle systems with the complicated spectra and strong 
interaction in a very natural and compact way. The essence of the method deeply 
related with the notion of the Generalized Mean Fields (GMF), which determine 
the elastic scattering corrections. The inelastic scattering corrections leads to the 
damping of the quasiparticles and are the main topic of the present consideration. 
The calculation of the damping has been done in a self-consistent way for both 
limits. For the random Hubbard model the weak coupling case has been considered 
and the self-energy operator has been calculated using the combination of the IGF 
method and Coherent Potential Approximation (CPA). The other applications of 
the method to s-f model, Anderson model, Heisenberg antiferromagnet, electron- 
phonon interaction models and quasiparticle tunneling are discussed briefly. 
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1 Introduction 



The study of the Highly Correlated Electron Systems has attracted much attention re- 
cently, especially after discovery of copper oxide superconductors and the new class of 
heavy fermion compounds, coexisting with magnetism |ll - 0. Although much work has 
been performed during last years it is worthy to remind that the investigation of the exci- 
tations in many-body systems has been one of the most important and interesting subject 
for last few decades. The quantum field theoretical techniques have been widely applied 
to statistical treatment of a large number of interacting particles. Many-body calculations 
are often done for model systems of statistical mechanics using perturbation expansion. 
The basic procedure in many-body theory |^ is to find relevant unperturbed Hamiltonian 
and then take into account the small perturbation operator. This procedure, which work 
well for the weakly interacting systems, needs the suitable reformulation for the many- 
body systems with complicated spectra and strong interaction. For many practically 
interesting cases the standard schemes of perturbation expansion must be reformulated 



greatly P] - [|T2|. The most characteristic feature of the recent advancement in basic 
research on electronic properties of solids is development of variety of the new class of 
materials with unusual properties: high superconductors, heavy fermion compounds, 
diluted magnetic semiconductors etc. Contrary to the simple metals, where the funda- 
mentals very well known and the electrons can be represented in a way such that they 
weakly interact with each other (c.f. in these materials the electrons interact strongly 
and moreover their spectra are complicated, i.e. have many branches etc. This gives rise 
to interesting phenomena such as magnetism, metal-insulator transition in oxides, heavy 
fermions etc., but the understanding of what is going on is in many cases only partial 
if exist at all. Therefore the theoretical studies of the Highly Correlated Electron Sys- 
tems (HCS) are very important and actual. A principle importance of of these studies is 
concerned with a fundamental problem of electronic solid state theory, namely with the 
tendency of 3d electrons in TMC and 4f electrons in rare-earth metal compounds (REC) 
and alloys to exhibit both localized and delocalized behaviour. The interesting electronic 
and magnetic properties of these substances are intimately related to this dual behaviour 
of electrons. In spite of experimental and theoretical achievements 111 - 0, still it remains 
much to be understood concerning such systems. Recent theoretical investigations of HCS 
have brought forth significant variety of the approaches which are trying to solve these 
controversial problems. It seems appropriate to point out that a number of perturbation- 
theory or mean-field theory approaches which have been proposed in the past few years, 
are in fact questionable or inadequate. In order to match such a trend we need to develop 
a systematic theory of the Highly Correlated Systems, to describe from the first principles 
of the condensed matter theory and statistical mechanics the physical properties of this 
class of the materials. In the present paper we will present the approach which allows one 
to describe completely the quasi-particle spectra with damping in a very natural way. This 
approach has been suggested as essential for various many-body systems and we believe 



that it bear the real physics of Highly Correlated Systems [|I^ , |T^ . The essence of our 
consideration of the dynamical properties of many-body system with strong interaction 
is related closely with the field theoretical approach and use the advantage of the Green's 
functions language and the Dyson equation. It is possible to say that our method tend to 
emphasize the fundamental and central role of the Dyson equation for the single-particle 
dynamics of the many-body systems at finite temperature. 



2 



Just this point differ our IGF approach from the comphmentary many-body approach 
which is based on the moment expansion for the spectral functions. It was developed in 
a very detail by W.Nolting - . 



2 Irreducible Green's Functions Method 

In this Section, we will discuss briefly the novel nonperturbative approach for description 
of the many-body dynamics of the HCES. At this point it is worthwhile to underline that 
it is essential to apply an adequate method in order to solve a concrete physical problem; 
the final solution should contain a correct physical reasoning in a most natural way. The 
list of many-body techniques that have been applied to strongly correlated systems is 
extensive. The problem of adequate description of many-body dynamics for the case of 
very strong Coulomb correlations has been explicitly raised by Anderson, who put the 
direct question: "... whether a real many-body theory would give answers radically dif- 
ferent from the Hartree-Fock results?" |^ (c.f. [^). The formulation of the Anderson 
modeipri and closely related Hubbard model PSI, [El] dates really a better understand- 



ing of the electronic correlations in solids, especially if the relevant electrons are modelled 
better by tight-binding approximation [^,||2B|. Both of the models, Anderson and Hub- 
bard, are often referred to as simplest models of magnetic metals and alloys. This naive 
perception contradicts the enormous amount of papers which has been publishing during 
the last decades and devoted to attacking the Anderson/Hubbard model by many refined 
theoretical techniques. As is well known now, the simplicity of the Anderson/Hubbard 
model manifest itself in the dynamics of a two-particle scattering. Nevertheless, as to 
the true many-body dynamics, there is still no simple and compact description. In this 
paper it will be attempted to justify the use of a novel Irreducible Green's Functions 
(IGF) for the interpolation solution of the single-band Hubbard model and other basic 
solid state models as s-f model [^, Anderson model Heisenberg antifer- 

romagnet [^] and strong electron-phonon interaction model in modified tight-binding 
approximatiom(MTBA) for normal and superconducting metals and alloys ||33|, |34 



A number of other approaches has been proposed and the our approach is in many re- 
spect an additional and incorporate the logic of development of the many-body techniques. 
The considerable progress in studying the spectra of elementary excitations and thermo- 
dynamic properties of many-body systems has been for most part due to the development 
of the temperature dependent Green's Functions methods. We have developed the help- 
ful reformulation of the two-time GFs method which is especially adjusted [^] for the 



correlated fermion systems on a lattice. The very important concept of the whole method 
are the Generalized Mean Fields. These GMFs have a complicated structure for the 
strongly correlated case and are not reduced to the functional of the mean densities of 
the electrons, when we calculate excitations spectra at finite temperature. To clarify the 
foregoing, let us consider the retarded GF of the form 

=« A{t), Bit') »= -lOit - f) < {A{t)B{t% >,r] = ±1. (1) 

As an introduction of the concept of IGFs let us describe the main ideas of this approach 
in a symbolic form. To calculate the retarded GF G{t — t') let us write down the equation 
of motion for it: 

uG{u}) =< [A, A+]r, > + « [A, H]_\A+»^. (2) 
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The essence of the method is as follows flM. It is based on the notion of the "IRRE 



DUCIBLE" parts of GFs (or the irreducible parts of the operators, out of which the GF 
is constructed) in term of which it is possible, without recourse to a truncation of the 
hierarchy of equations for the GFs, to write down the exact Dyson equation and to obtain 
an exact analytical representation for the self-energy operator. By definition we introduce 
the irreducible part (ir) of the GF 

« [A, H]_\A+ »=« [A, H]_ - zA\A+ » . (3) 

The unknown constant z is defined by the condition (or constraint) 

<[[A,Ht,A\>=0 (4) 

From the condition (4) one can find: 

^ < [[A,H].,A+],> ^ Ml 
' <[A,A+],> Mo 

Here Mq and Mi are the zeroth and first order moments of the spectral density. Therefore, 
irreducible GF (3) are defined so that it cannot be reduced to the lower-order ones by any 
kind of decoupling. It is worthy to note that the irreducible correlation functions are well 
known in statistical mechanics. In the diagrammatic approach the irreducible vertices are 
defined as the graphs that do not contain inner parts connected by the G°-line. With the 
aid of the definition (3) these concepts are translating into the language of retarded and 
advanced GFs. This procedure extract all relevant (for the problem under consideration) 
mean field contributions and puts them into the generalized mean-field GF, which here 
are defined as 

G%cu) = ^ ^^^^ ^ . (6) 

To calculate the IGF << [A, H]_{t), A'^{t') » in (2), we have to write the equation 
of motion after differentiation with respect to the second time variable t'. The condition 
(4) remove the inhomogeneous term from this equation and is the very crucial point of 
the whole approach. If one introduces an irreducible part for the right-hand side operator 
as discussed above for the "left" operator, the equation of motion (2) can be exactly 
rewritten in the following form 

G = G° + G°PG°. (7) 

The scattering operator P is given by 

p = (Mo)"^ « [A, H]_\[A+ , H]_ {Mo)-\ (8) 

The structure of the equation (7) enables us to determine the self-energy operator M, in 
complete analogy with the diagram technique 

P = M + MG^P. (9) 

From the definition (9) it follows that we can say that the self-energy operator M is defined 
as a proper (in diagrammatic language "connected" ) part of the scattering operator M = 
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{Py. As a result, we obtain the exact Dyson equation for the thermodynamic two-time 
Green's Functions: 

G = G^ + G^MG, (10) 
which has well known formal solution of the form 

G= [{G^)-^-M]-\ 

Thus, by introducing irreducible parts of GF (or the irreducible parts of the operators, out 
of which the GF is constructed) the equation of motion (2) for the GF can be exactly (but 
using constraint (4)) transformed into Dyson equation for the two-time thermal GF. This 
is very remarkable result, which deserve the underlining, because of the traditional form of 
the GF method did not included namely this point. The projection operator technique 



has essentially the same philosophy, but with using the constraint (4) in our approach we 
emphasize the fundamental and central role of the Dyson equation for the calculation of 
the single-particle properties of the many-body systems. It is important to note, that for 



the retarded and advanced GFs the notion of the proper part is symbolic in nature ||14 
However, because of the identical form of the equations for the GFs for all three types 
(advanced, retarded and causal), we can convert in each stage of calculations to causal GFs 
and, thereby, confirm the substantiated nature of definition (9)! We therefore should speak 
of an analog of the Dyson equation. Hereafter we will drop this stipulation, since it will not 
cause any misunderstanding. It should be emphasized that scheme presented above give 
just an general idea of the IGF method. The specific method of introducing IGFs depends 
on the form of operator A, the type of the Hamiltonian and the conditions of the problem. 
The general philosophy of the IGF method lies in the separation and identification of 
elastic scattering effects and inelastic ones. This last point is quite often underestimated 
and both effects are mixed. However, as far as the right definition of quasiparticle damping 
is concerned, the separation of elastic and inelastic scattering processes is believed to be 
crucially important for the many-body systems with complicated spectrum and strong 



interaction. Recent paper [0 emphasizes especially that the anomalous damping of 



electrons (or holes) distinguishes cuprate superconductors from ordinary metals. From a 
technical point of view the elastic (GMF) renormalizations can exhibit a quite non-trivial 
structure. To obtain this structure correctly, one must construct the full GF from the 
complete algebra of the relevant operators and develop a special projection procedure 
for higher-order GF in accordance with a given algebra. The Hubbard model is a very 
suitable tool for the applying of this approach 



3 Hubbard Model 

The model Hamiltonian which is usually referred to as Hubbard Hamiltonian 

H = ^UjafaO-ja + U/2j2'n'iani-a (11) 

ij(T id 

includes the intraatomic Coulomb repulsion U and the one-electron hopping energy tij. 
The electron correlation forces electrons to localize in the atomic orbitals, which are 
modelled here by the complete and orthogonal set of the Wannier wave functions [0(r — 
Rj)]- On the other hand, the kinetic energy is reduced when electrons are delocalized. The 
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main difficulty of the right solution of the Hubbard model is the necessity to taking into 
account of the both these effects simultaneously. Thus, the Hamiltonian (11) is specified 
by two parameter: U and effective electron bandwidth 

ij 

The band energy of Bloch electrons e{k) is defined as follows 

k 

where the N is the number of the lattice sites. It is convenient to count the energy from 
the center of gravity of the band, i.e. = J2k ^ (^) = 0- The effective electron bandwidth 
A and Coulomb intrasite integral U define completely the different regimes in 3 dimension 
depending on parameter 7 = A/f/. It is usually a rather difficult task to find interpolation 
solution for the dynamical properties of the Hubbard model. To solve this problem with 
a reasonably accuracy and describe correctly an interpolating solution from "band" limit 
(7 S> 1) to "atomic" limit (7 — >• 0) one need more sophisticated approach than usual 
procedures which have been developed for description of the interacting electron-gas- 
problem. We evidently have to to improve the early Hubbard's theory taking account 
of variety of possible regimes for the model depending on electronic density, temperature 
and values of 7. The single-electron GF 

Gij^iuj) =« a^^\a^^ »= A^"^^G<,(A;,u;)exp[-iA;(^i - Rj)], (12) 



which has been calculated by Hubbard ||23|, has the characteristic two-pole functional 
structure 

G^{k,cu) = [F^ico)-eik)]-' (13) 

where 

" ^ ' {uj - E+ - nZ^X){u - E^ - nt^X) - nt^nZ^X^ ^ ' 

andA is the certain function which depends on parameters of the Hamiltonian. If A is 
small (A —>■ 0) then expression (14) take the form: 

uo-E^- nt^X ^ uj-E+-nZ,X' 
which correspond to the two shifted subbands with the gap 

uJi-uo2 = iE+ - E_) + {nZ^ - n+JA = U + A2n+^. 
Here n"*" = n and n~ = 1 — n; E^ = U, E_ = 0. If A is very big then we obtain 

A 1 



Fzh 



UJ 



[{io - E^)n_„ + E+)nt^]X uj - {nt^E+ - n_^E. 



This latter solution correspond to the single band, centered at the energy uo n'tJJ. 
The two- pole functional structure of the single-particle GF is very easy to understand 
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within formalism which describe the motion of electrons in binary alloys [40|. If one 
introduce the two types of the scattering potentials t± ^ {u — E±)~^ then the two kinds 
of the t-matrix and T_ appears which satisfy the following system of equations: 



where is the bare propagator between the sites with the energies E±. The solution 
of this system has the following form 



t± + t±G^t± 



(1 - t+G\+){l - t-G'L.) - G\Gl_t+t- 

t^' + Gl 

(t;' - Gl+)itZ^ - GO _) - G\Gl. 



(15) 



Thus, by comparing this functional two-pole structure and well-known "Hubbard III" 
solution 



Tj^{uj) = uj - F^{uj) 

it is possible to identify the "scattering corrections" and "resonance broadening correc- 
tions" in the following way: 



uj{uj -U) -{uo - Un_^)A^{ijL!) 



k 



[UJ 



If we put Acr{Lu) = we immediately obtain the "Hubbard I" solution p3|]. The "al- 
loy analogy" approximation correspond to Aa-^u) ^ Y„{uj). Note, that the "Hubbard 
IH" self-energy operator Ti^{uj) is local, i.e. do not depend on quasimomentum. The 
another drawback of this solution is very inconvenient functional representation of the 
elastic and inelastic scattering processes. The conceptually new approach to the theory 
of very strong but finite electron correlation for Hubbard model has been proposed by 



Roth 



She clarified microscopically the origination of the two-pole solution of the 
single-particle GF, what was the very unusual fact from the point of view of the standard 
Fermi-liquid approach, showing that the naive one-electron approximation of the band 
structure calculations is not valid for the description of the electron correlations in HCS . 
Thus the use of sophisticated many-body technique is required for the calculation of the 
excitation spectra at finite temperature. This last point should be underlined, because 



of the suitable modification of the Density Functional Approximation could 
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give the reasonable description of the ground-state properties of HCS. We shall show 
here, following the papers ^^at the use of the IGF method permit to improve 

substantially both solutions, Hubbard's and Roth's, by defining the correct Generalized 
Mean Fields for the Hubbard model. 



4 Hubbard Model. Weak Correlation 

The concept of the GMFs and the relevant algebra of operators from which GFs are 
constructed are the central ones to our treatment of electron correlation in solids. It will 
be convenient (and much more shorter) to discuss these concepts for weakly and strongly 
correlated cases separately. For the first time we must to construct the suitable state 
vector space of the many-body system The fundamental assumption implies that 

the states of a system of interacting particles can be expanded in terms of the states of 
non-interacting particles This concept originate in perturbation theory and finds 



support for weakly interacting many-particle systems(c.f. 0). For the strongly correlated 
case this approach needs the the suitable reformulation (cf. [Q) and namely in this point 
the right definition of the GMFs is vital. Let us consider the weakly correlated Hubbard 
model (11). In many respect this case is similar to the ordinary interacting electron gas 
but with very local, singular interaction. It will be shown below that the usual creation 
af^ and annihilation second quantized operators with the properties 

a+vl>(0) = vl/«;a,vl/(i) = vl>(0) 

a,vl/(0)=0;a,vl/f) = 0, ^ j) 

are suitable variables for the description of the considering systems. Here \1/ and '^^^^ are 
the vacuum and single-particle states respectively. The question now is how to describe our 
system in terms of the quasiparticles. For a translationally invariant system, to describe 
the low-lying excitations of the system in terms of quasiparticles one has to choose 



eigenstates such that they all correspond to definite momentum. For the single-band 
Hubbard model (11) the exact transformation reads 

i 

Note, that for degenerate bands model the more general transformation is necessary. Then 
the Hubbard Hamiltonian (11) in the Bloch vector state space are given by 

H = Y,^{k)at„aka + t//2A^ ^ ^ a+.^_g„ap^a+_^a^_^ (16) 

If the interaction is weak, the algebra of the relevant operators is very simple: it is an 
algebra of the non-interacting fermion system (0^0-, cil^, n^^ = a'^^a^a)- For the calculation 
of the electronic quasiparticle spectrum of the Hubbard model in this limit let us consider 
the single-electron GF, which are defined as 

Gkait - f) =« aka, at »= - t') < [aka{t), a+^{t')]+ >= 

/+00 
du exp{—iujt)Gkai^^) = 
-00 

r+co r+co d^j' 

1/271 dujexp{-iujt)l/27T -{exp{f3uj') + l)Akaiuj') (17) 

J —00 J— 00 UJ — UJ 



8 



where (5 — {kT)^i and Ak„{u!) is the spectral intensity. The equation of motion for the 
Fourier transform of the GF Gkcr{u;) has the form 

{uj - ek)Gka{'^) = 1 + U/NY^ « ak+paa^+g_^ag-a\at^ >>a; (18) 

pq 

Let us introduce, by definition, an "irreducible" GF in the following way 

<< (^k+paO-p+q-cr^^q-alo-ta ~^p,0 < f^-q-a > Gka (19) 

The irreducible (ir) GF in (19) is defined in such a way that it cannot be reduced to GF 
of lower order with respect to the number of fermion operators by an arbitrary pairing of 
operators or, in another words, by any kind of decoupling. Substituting (19) in (18) we 
obtain 

Gfe,(cu) = G^f (a;) + G^^(^)[//Ar^ - « afc+p,a++,_,a,_,|a+ (20) 

pq 

Here we have introduced the notations 

Gff{uj) ^{u- e{ka))-'; e{ka) = e{k) + U/NY.< > (21) 

In this paper, for brevity, we confine ourself by considering the paramagnetic solutions 
only, i.e. < >=< ri-a >■ In order to calculate the higher-order GF on the r.h.s. of 
(20) we have to write the equation of motion obtained by means of differentiation with 
respect to the second variable f. Constraint (4) allows us to remove the inhomogeneous 
term in this equation for -^^^ « A{t), a^^(t') >>. 
For the Fourier components, this is written in the form 



{u - e{k)y' « A\at »a;=<'' [A, a+ ]+ > + 

<< A\a'^_^aj.+s-a(J'k+sa »ui ■ (22) 

rs 

The anticommutator in (22) is calculated on the basis of the definition of the irreducible 
part 

< [«fc+pCTap+q-(T«q-(T- < dp+q-af^q-a > 0-k+pa, (^1^]+ ^ (23) 

If one introduces irreducible part for the r.h.s. operators by analogy with expression (19), 
the equation of motion (20) can be exactly rewritten in the form (7) 

GU^) = G^fiuj) + Grf{u;)PU^)Gr/{uj) (24) 

where we have introduced the following notation for the operator P (8) 

pqrs 

TP 

~ "a72 2^ "^"^ 0'k+p<j0'p+q^cr^q-a\0'r-cr(^r+s-<j0'k+s<7 »u) l^^J 
pqrs 
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To define the self-energy operator according to tlie (9) one must separate tlie "proper" 
part by tlie following way 

DLip, q\r, s; uj) = L]^^{p, q\r, s; u) 
+^ E LUp,q\r's';u;)G'j^^{u;)DUp',q'\r,s-,u;) (26) 

r's'p'q' 

Here L^^(p, g|r, s; u;) is the "proper" part of the GF D^^^{p,q\r, s]uj), which in accor- 
dance with the definition (19) cannot be reduced to the lower-order one by any type of 
decoupling. Using (9) we find 

Gfc. = GZ^iu) + Gti''{u)M^{k,uj)G,Au) (27) 

Equation (27) is the Dyson equation for the single-particle two-time thermal GF. Accord- 
ing to (10) it has the formal solution 

GkA^) = [lo- e{ka) - M,(fc, a)]-' (28) 

where the self-energy operator M is given by 



2 



Ma{k, w) = E ^tiP, q\r, s; u) 



iV2 

pqrs 

<< 0'k+pcrO'p+q-crO'q~a\0'k+scr(^r-cr(^r+s-cr » ) 



pqrs 



The latter expression (29) is an exact representation (no decoupling has been made till 
now) for the self-energy in terms of higher-order GFs up to second order in U (for the 
consideration of the higher order equations of motion see Ref. |]15|). Thus, in contrast 
to the standard equation-of-motion approach the determination of the full GF has been 
reduced to the calculation of the mean-field GF G*^^ and the self-energy operator M. 
The main reason for this method of calculation is that the decoupling is only introduced 
into self-energy operator, as it will shown in a detail below. The formal solution of the 
Dyson equation (28) define the right reference frame for the formation of the quasipar- 
ticle spectrum due to the its own (formal solution) correct functional structure. In the 
standard equation-of-motion approach such a structure could be lost by using decoupling 
approximations before arriving to the correct functional structure of the formal solution 
of the Dyson equation. This is a crucial point of the IGF method. The energies of the 
electronic states in the mean-field approximation are given by the poles of G^^ (21)- Now 
let us consider the damping effects and finite lifetimes. To find an explicit expression for 
self-energy M (29), we have to evaluate approximately the higher-order GF in (21). It 
will be shown below that the IGF method can be used to derive the damping in a self- 
consistent way simply and more generally than other formulations. First, it is convenient 
to write down the GF in (29) in terms of correlation functions by using the well-known 



spectral theorem [^] : 



O'k+pa0'p^q-.cr^q-o-\0'k+sa'^r-cr'-''r+s-cr 



/ + CXD 
exp(icj't) 
-oo 



2tt J-oo uj — uj' 

< al^^^{t)a'^_„{t)ar+s-a{t)ak+paap^q-aaq-a > (30) 
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Further insight is gained if we select the suitable relevant "trial" approximation for the 
correlation function on the r.h.s. of (30). In this paper we show that the earlier for- 
mulations, based on the decoupling or/and diagrammatic methods can be arrive at from 
our technique but in a self- consistent way. Clearly that the choice of the relevant trial 
approximation for correlation function in (30) can be done in many ways. For example, 
the reasonable and workable one may be the following "pair approximation", which is 
especially good for the low density of the quasiparticles: 

< '^fc+scr(^)'^^-cr('^)'^r+s-o-(^)a.fc+pcrCtp+g_o-S-o- >*'"~ 

^k+s,k+p^r,q^r+s,p+q (31) 

Using (30) and (31) in (29) we obtain the approximate expression for the self-energy op- 
erator in a self-consistent form (the self-consistency means that we express approximately 
the self-energy operator in terms of the initial GF and, in principle, one can obtain the 
required solution by suitable iteration procedure): 



UJ + UJi — UJ2 — ^3 



pq •' 

[n{uj2)n{uj3) + n{uJi){l - n{uj2) - n{uj3))]gp+g^^{uji)gk+p^{uj2)gq-a{^^3) (32) 

where we have used the notations 

1 



gkai^) = ImGkaiuj + ie);n{iu) = [exp{/3uj) + 1] 



-1 



77 



The equations (28) and (32) form a closed self-consistent system of equations for the single- 
electron GF for the Hubbard model, but for weakly correlated limit only. In principle, we 
may use, on the r.h.s. of (32) any workable first iteration-step form of the GF and find a 
solution by repeated iteration. It is most convenient to choose as the first iteration step 
the following simple one-pole approximation: 

gk.{uj) ^ 5{uj - e{ka)) (33) 

Then, using (33) in (32), we get for the self-energy an explicit and simple expression 

T\/f (u \ ^ ^p+<7-cr(l '^fc+pcr fT'q-cr) "I" ^k+pcr'^q-CF ( QA^ 

MAt^,u;)-—X^ ^ + eip + qa) - e{k + pa) - e{qa) ^^^^ 

The numerical calculations of the typical behaviour of real and imaginary parts of the 
self-energy (34) have been performed for the model density of states of the FCC 
lattice. These calculations and many other - prove that the conventional one- 
electron approximation of the band theory is not always a sufficiently good approximation 
for transition metals like nickel. The simple formula (32) derived above for the self- 
energy operator are typical in showing the role of correlation effects in the formation of 
quasiparticle spectrum of the Hubbard model. It is instructive to examine other types of 
the possible trial solutions for the six-operator correlation function in the eqn.(30). The 
approximation which we propose now reflects the interference between the one-particle 
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branch of the spectrum and the collective one: 



< a+_^(t)ag_^ >< a^+^^(t)ar+s-a{t)ak+paap+q-cr > (35) 

It is visible now that the three contributions in this trial solution describe the self-energy 
corrections that take into account the collective motions of electron density, the spin den- 
sity and the density of "doubles" , respectively. The essential feature of this approxima- 
tion is connected with the fact that correct calculation of the single-electron quasiparticle 
spectra with damping require the suitable incorporating of the influence of the collective 
degrees of freedom on the single-particle ones. The most interesting contribution is re- 
lated with the spin degrees of freedom because of correlated system are the magnetic or 
have very well developed magnetic fluctuations. We follows the above steps and calculate 
the self-energy operator (29) as 

M^{k,uj) = — dujidu2- ^ ^ ^ 



N J-oD U — UOi — U2 

^exp[-ik{Ri - Rj)]{--Im « S^\SJ 



TT 

where the following notations have been used: 

N{uj) = [exp{(3uj) - 
It is possible to rewrite (37) in a more convenient way now 



^ Im « ai_^\aj_^ >>i^2) (36) 



M„{k,u) = —Yl J ^^(cot^^ + tan— ^ 



{--Imx'^^{k-q,u-uj')gq^{uj')). (37) 

TT 

The equations (28) and (37) form again another self-consistent system of equations for the 
single-particle GF of the Hubbard model. Note, that both expressions for the self-energy 
depend on quasimomentum; in other words the approximate procedure do not broke the 
momentum conservation law. It is important, because of the poles uj{k, a) = e{k, a)—ir{k) 
of the GF (28) are determined by the equation 

u - e{ka) - Re[M„{k, cu)] - iIm[M^{k, cu)] = (38) 

It may be shown quite generally that the Luttinger's definition of the true Fermi surface 
is valid in the framework of our theory. It is worthy to note that for electrons in a crystal 
where there is a band index, as well as quasimomentum, the definition of the Fermi surface 
are a little more complicated then the single-band one. Before the single particle energies 
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and Fermi surface are known, one must carry out a diagonalization in the band index. In 
order to give a complete picture of the GMFs let us discuss briefly the interesting question 
of the correct deflnition of the so-called unrestricted Hartree-Fock approximation (UHFA). 
Recently, this approximation has been applied for the single-band Hubbard model (11) 
for the calculation of the density of states for Cu02 clusters [0. The following deflnition 
of UHFA has been used: 

Thus, in addition to the standard HF term, the new, the so-called "spin-flip" terms, 
are retained. This example clearly show that the nature of the mean-flelds follows from 
the essence of the problem and should be deflned in a proper way. It is clear, however, 
that the deflnition (39) broke the rotational symmetry of the Hubbard Hamiltonian. For 
the single- band Hubbard Hamiltonian the averaging < af_^ai^a >= because of the 



rotational symmetry of the Hubbard model. So, in Ref. |51[ the effective Hamiltonian 
Hes has been deflned. We have analysed in detail the proper deflnition of the irreducible 
GFs which include the "spin-flip" terms. The deflnition (19) must be modifled in the 
following way: 

5pfl < Uq^^ > Gka- < ak+paap^q_^ >« a^-^ja^^ >>^ (40) 

From this deflnition follows that such a type of introduction of the IGF broaden the initial 
algebra of the operator and initial set of the GFs. That means that "actual" algebra of 
the operators must include the spin-flip terms at the beginning, namely: (ojo-, af^, rii^, 
af^ai^^). The corresponding initial GF will have the form 

<< ai^\a^^ » « ai„\aj_^ » 
« ai_^\aj^ » « ai_^\aj_^ » 

In fact, this approximation has been investigated earlier by Kishore and Joshi ||52 
They clearly pointed out that they assumed that the system is magnetized in x direction 
instead of conventional z axis. 



5 Hubbard Model. Strong Correlation 

When studying the electronic quasiparticle spectrum of the strongly correlated systems, 
one must take care of at least three facts of major importance: 

(i) The ground state is reconstructed radically as compared with the weakly correlated 
case. Namely this fact lead to the necessity of the redeflnition of the single-particle states. 
Due to the strong correlation, the initial algebra of the operators are transformed into 
new algebra of the complicated operators. In principle, in terms of the new operators the 
initial Hamiltonian may be rewritten as bilinear form and the generalized Wick theorem 



can be formulated [Q, [Q]. It is very important to underline, that the transformation to 
the new algebra of relevant operators reflect some important internal symmetries of the 
problem and nowadays this way of thinking are formulating in elegant and very powerful 



technique of the classiflcation of the integrable models , ||56[ and exactly soluble models 
(see also|p7|). 
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(ii) The single-electron GF, which describe the dynamical properties, must have two-pole 
functional structure, giving in the atomic limit, when hopping integral tends to zero, the 
exact two-level atomic solution. 

(iii) The GMFs have, in general case, a very non-trivial structure. The GMFs functional 
cannot be expressed in terms of the functional of the mean particles density. 

In this section we consider large, but finite. Coulomb repulsion. The inspiring ideas of 
papers |3^, where the problem of the relevant algebra of the operators has been 

considered, are central to our consideration here. Following this approach we consider the 
new set of relevant operators: 



nt^a,^, (a = ±); n+ = n^^, n^^ 



'1 - UiJ: 



(41) 



a 

The new operators diaa and (i^^^ have complicated commutation rules, namely 

The convenience of the new operators follows immediately if one write down the equation 
of motion for them 



(42) 



It is possible to interpret both contribution in this equation as alloy analogy 

and resonance broadening correction. Let us consider the single-particle GF (12) in the 
Wannier basis. Using the new operator algebra it is possible to rewrite identically GF 
(12) in the following way 



af3 



The equation of motion for the auxiliary matrix GF 



have the following form 

(EFij„{uj) - I6ij)ap = « n1_^ai^ + aai„hii-„\d'^p^ >>, 

Where the following matrix notations have been used 



E 



[uj - E+) 
iuj-E^ 



I 



n 



nl 



(43) 



(44) 



(45) 



(46) 



In accordance with the general method of Section 2 we introduce by definition the matrix 
IGF: 



A+a' 
A -a' 



L IIO" iio" J 



(47) 



14 



Here the notations have been used: 

It is worth to underhne that the definition (47) are in heart of the whole our approach to 
description of the strong correlation in the Hubbard model. The coefficients A and B are 
determined from the constraint (4), namely 

<[(Di[^.)«/3,4/3j+>=0 (48) 
After some algebra we obtain from (48) {i ^ j) 

[Ai]ai3 = a(< djp_^ai-^ > + < di_p_^al_^ >){nt^)~^ 

< ai^ai-„al_^al^ >)]{n^^y^ (49) 

As previously, we introduce now the GMF GF F?^-^ in analogy with (6), however, as it is 
clear from (47), the actual definition of the GMF GF is very non-trivial. After the Fourier 
transformation we get 



Fl-)- ah-cd\nt^c nZ^a ^ 
The coefficients a, 6, c, d are equal to 

p 

':=N-'Y: e{p){A^^{-p) - B^^ip - q)) (51) 
a p 

Then, using the definition (43) we find the final expression for the GMF GF 

G^^F^k uj) = ^ - jnt^E. + nZ^E+) - \{k) ^^^^ 

{u - E+ - nZ„Xi{k)){i^ - - nt„\2{k)) - nZ„nt„\^{k)Xi{k) 



Here we have introduced the following notations: 
\,{k) 1 



Xsik) 1 



A4(fc) n"^ 



Y^e{p){A^^{-p)-B^^{p-k)) (54) 



A(A:) = (n:j2(Ai + A3) + (nt^nX^ + A4) 

From the equation (52) it is obvious that our two-pole solution is more general than "Hub- 
bard HI" 1^ and Roth[^ solutions. Our solution has the correct non-local structure. 



taking into account the non-diagonal scattering matrix elements more accurately. Those 
matrix elements describe the virtual "recombination" processes and refiect the extremely 
complicated structure of the single-particle state, which virtually include a great number 
of intermediate scattering processes (c.f. interesting analysis in Ref. ||55|). 
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The spectrum of the mean-field quasiparticle excitations follows from the poles of the GF 
(52) and consist of two branches 



uj^{k) = 1/2[{E+ -E^+ai + hi) ± ^{E+ + E^ - - hif - Acd] (55) 

where ai(6) = uj — E± — a{h). Thus the Spectral Intensity Aka{oj) of the GF (52) consist 
of two peaks, which separated by the distance 

ai — hi 



^1 - ^2 = V(f/ - ai - hiY -cd^ U{1 - ^-j^) + 0(7) (56) 

For the deeper insight into the functional structure of the solution (52) and to compare 
with the other solutions we rewrite the (50) in the following form 




from which we obtain for the {k^ u 



n 



+ (l + c6-i) nzJl + da-^] 



Gr{k,u) = -r-z.' + 



a — dh h — ca ^d 



u - E_-nl,Wl„{k) UJ - E+ - nZ„W+„{k) 



(58) 



where 



((< atantaj-<T > + < ai^„nf^a+_^ >) + 
(< nf^^nf^^ > + < ai^af_^aj-aafa > - < ai^ai_„aj_^aj^ >)) (59) 

are the shifts for the upper and lower splitted subbands due to the elastic scattering of 
the carriers in the Generalized Mean Field. Namely are the functionals of the GMF. 
The most important feature of the present solution of the strongly correlated Hubbard 
model is a very nontrivial structure of the mean- field renormalizations (59), which is 
crucial to understanding the physics of strongly correlated systems. It is important to 
emphasize that namely this complicated form of the GMF are only relevant to the essence 
of the physics under consideration. The attempts to reduce the functional of the GMF 
to the simpler functional of the average density of electrons are incorrect namely from 
the point of view of the real nature of the physics of HCS. This physics clearly show that 
the mean-field renormalizations cannot be expressed as a functional of the electron mean 



density. To explain this statement let us derive the "Hubbard 1" solution I^S] from our 
GMF solution (52). If we approximate (59) as 



nt^nZ^W^ik) ^ N-'J2Ujexp{-ik{Ri - Rj)) < nf_^nt^ > (60) 



and makes the additional approximation, namely 
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then solution (52) goes over into the "Hubbard I" solution 

^) ^ + )^ ^ (61) 

uj — U — e{k)n_a- OJ — e{k)(l — n_cr) 

This solution, as it is well known, is unrealistic from the many points of view. 

As regards to our solution (52), the second important aspect is that the parameters 

\i{k) do not depend on frequency, i.e. depends essentially on the elastic scattering pro- 
cesses. Such a dependence on frequency arises due to inelastic scattering processes which 
are contained in our self-energy operator and we proceed now with the derivation of the 
explicit expression for it. 

To calculate the high-order GF on the r.h.s. of (45) we should use the second time variable 

{t') differentiation of it again. If one introduces irreducible parts for the right-hand-side 
operators by analogy with expression (47), the equation of motion (45) can be rewritten 
exactly in the following form 

= + (62) 

Here the scattering operator P (8) has the form 

P,.((^) = l-'iY^tat^, « D:[JD-+^. (63) 

Im 

In accordance with the definition (9) we write down the Dyson equation 

F = F° + F°MF (64) 

The self-energy operator M is defined by Eq. (9). Let us note again that the self-energy 
corrections, according to (10), contribute to the full GF as an additional terms. This is 
an essential advantage in comparison with the "Hubbard III" solution and other two-pole 
solutions. For the full GF we find, using the formal solution of Dyson equation 

+^^d + Ti^M'iK 00)) + ^{c + nZ^M;+{k, u))\ 



[det((F°,(u;))-^-M.(^,u;))]-^ (65) 



After some algebra we can rewrite this expression in the following form, which is essentially 
new and, in a certain sense, are the central result of the present theory 

Q ^ uJ-{n+E_+n-E+)-L 

(oj — — n~Li)(u; — — n+L2) — n'n+LsL^ 



where 

L,{k,u) = X,{k) - ■^M+'-{k,u); 
L^iKu) = X2{k) - ^M--{k,u); 
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T) 



Tl 

Li{k,u) = A4(A;) + — ^M-+(A;,a;); 



n_ ^ 

L{k, oj) = X{k) + nl^n:^(M++ + M~ - M-+ - M+-) (67) 

Thus, now we have to find the exphcit expressions for the elements of the self-energy 
matrix M. To proceed we should use the spectral theorem again, as in Eq. (30), to 
express the GF in terms of correlation functions 

M:'f^{k,uj)^<D^:^l^{t)Dl^> (68) 

For the approximate calculation of the self-energy we propose to use the following trial 
solution 

< L>-+(i)L>- >^< ai,{t)ai, X nl,{t)r^_, > 
+ < a+„{t)nf_, >< n^_^(t)ai^ > +(3 < bij_^(t)ai^ >< a%{t)nt^ > 

+/3 < b+j_^{t)n'l_^ >< aj„{t)ai„ > +« < a+^{t)ai^ >< n!^_„it)hu^a > 

+a < a:^„{t)bii_^ >< n^_^{t)bu.^ > 
+a(3 < b^j_„(t)aia >< a^„{t)bii-„ > 

< bi^_^{t)bii_, >< a+^{t)ai„ > (69) 

It is quite natural to interpret the contributions in this expression in terms of scattering, 
resonance-broadening and interference corrections of different types. For example, let us 
consider the simplest approximation. For this aim we retain the first contribution in (69) 

/+°° fid/ 
-^{eMPco') + 1) 
-oo 00 — UJ 



+00 

— exp{iuj't)N~^ V exp{-ik{Ri - Rj))tiitmj 



L 

J dujin{uji) exp{iuit)gmia{i^i) (^~^^^?f i'^i - • (70) 

Equations (70) and (64) are the self-consistent system of equations for the single-particle 
Green's function. For a simple estimation, for the calculation of the self-energy (70) it 
is possible to use any initial relevant approximation of the two-pole structure. As an 
example we take the expression (61). We then obtain 

[IMlU«^|6(fc-g)pi^f 
g 

!W . , . 

^uj-U-e{k-q)n-a uj-e{k-q){l-n-a) 

On the basis of the self-energy operator (71) we can explicitly find the energy shift and 
damping due to inelastic scattering of the quasiparticles, which is a great advantage of 
the present approach. It is clear from the present consideration that for the systematic 
construction of the approximate solutions we need to calculate the collective correlation 
functions of the electron density and spin density and the density of doubles, but this 
problem must be considered separately. 
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6 Correlations in Random Hubbard Model 



In this chapter we shall apply IGF method for consideration of the electron-electron 
correlations in the presence of disorder to demonstrate the advantage of our approach. 
The treatment of the electron motion in substitutionally disordered disordered A^Bi^x 
transition metal alloys is based upon certain generalization of Hubbard model, including 
random diagonal and off-diagonal elements caused by substitutional disorder in the binary 
alloy. The electron-electron interaction play an important role for various aspects of 
behaviour in alloys, e.g. in the weak localisation in Ti-Al alloys |]59[| (for recent review 
see [§n|)- There are certain aspects of the High-Tc superconductivity where disorder 
play a role and recently it have been discussed in papers^l[|, [^, where the distribution 
of magnetic molecular fields has been treated within the single-site Coherent Potrntial 
Approximation(CPA) The CPA has been refined and developed in many papers (e.g. 
|6^ , |65| ) and till now|66[| are the most popular approximation for theoretical studying of 



alloys. But the simultaneous effect of disorder and electron- electron inelastic scattering 
has been considered for some limited cases only and not within the self-consistent 



scheme. Let us consider the Hubbard model Hamiltonian on a given configuration of alloy 

^(.) _ ^(-) + ^(-) (72) 



where 

ia ija 

H^"^ = \Y.U>z<rni.„ (73) 

ia 

Contrary to the periodic model (11), the atomic level energy e^, the hopping integrals t"^ 
as well as the intraatomic Coulomb repulsion here are the random variables, which 
take the values , f^^ and W^, respectively; the superscript i/(/i) refers to the atomic 
species (i^, = A,B) located on site i(j). The nearest-neighbour hopping integrals are 
included only. 

To unify the IGF method and CPA into completely self-consistent scheme let us consider 
the single-electron GF (17) Gija in the Wannier representation for a given configuration 
[i/). The corresponding equation of motion has the form(for brevity we shall omit the 
superscript {p) where its presence is clear) 

[u - Ei) « ai„\aj^ »uj= Sij + ^tin « a^alaja »uj 

n 

+ Ui « rii-adialaja (^4) 

In the present paper, for brevity, we will confine ourselves by the weak correlation case 
and the diagonal disorder only. The generalization for the case of strong correlation or 
off-diagonal disorder is straightforward, but its length considerations preclude us from 
discussing at this time. 

Using the definition (3), we define the IGF for a given (fixed) configuration of atoms in 
an alloy as follows 



<< ni-aaia\a^^ »=« ni-aaia\a^^ » - < ni_„ >« ai<^|a+^ >> (75) 
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This time, contrary to (19), because of lack of translational invariance we must take into 
account the site dependence of < nj_o- >. Then we rewrite the equation of motion (76)in 
the following form 

Y^[{UJ -Ei-Ui < rii^^ >)6ij - tin] « ana\aj„ »uj = 
n 

6ij + Ui{"'« ni_^ai^\aj^ »^) (76) 

In accordance with the general method of Section 2, we find then the Dyson equation for 
a given configuration (z/) 

G.j^iu) = G° ,(.;) + Y: G',^,{u)Mmna{u;)Gnj.{u;) (77) 

mn 

The GMF GF G^-^ and the self-energy operator M are defined as 

"Y HimaG'l^j^{uj) = 5ij 
m 

Pmna -^mna ~(" ^ ^ ^^miaG j^j^Pjncr 

Hima (^ Ui <Z ^i—a ^)'^im ^im 

Pmnai'^) = « nm~aClmcT\nn~<T<4ia »Z)Un (78) 

In order to calculate the self-energy operator M self-consistently we have to express it 
approximately by the lower-order CPs. Employing the same pair approximation as (31) 
(now in Wannier representation) and the same procedure of calculations we arrive at the 
following expression for M for a given configuration (z/) 

ImG^nl-a {^i ) ImG^rl-a (^2 ) ImG^l^l^ (uj^); 
^ _ duJiduj2du3 (1 - n{uJi))n{u2)n{uj3) ^r^^^ 

As we have mentioned previously, all the calculations just presented have been done for 
a given configuration of atoms in alloy. All the quantities in our theory (G, G^, P, M) 
depends on the whole configuration of the alloy. To obtain a theory of a real macroscopic 
sample, we have to average over various configurations of atoms in the sample. The 
configurational averaging cannot be exactly made for a macroscopic sample. Hence we 
must resort to an additional approximation. It is obvious that self-energy M is in turn 
the functional of G, namely M = M[G]. If the process of taking configurational averaging 
is denoted by G, than we have 

G = G^ + CFMG 

Pew words are now appropriate for the description of general possibilities. The calculations 
of G° can be performed with the help of any relevant available scheme. In the present 
work, for the sake of simplicity, we choose the single-site CPApH], namely we take 



(, ,\ _ /\r-i exp(zfc(-Rm -Rn.)) 
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Here e{k) — in,o exp(i/ci?„), z is the number of nearest neigbours of the site 0, and 
the Coherent potential is the solution of the CPA self-consistency equations. For 

the AxBi_x these read 

E'^H = xe'\ + (1 - x)e% - (s^ - E'^)F-(a;, E-)(£^ - E'^); 

F'^(u;,E-) = G'^_H (81) 

Now, let us return to the calculation of the configurationally averaged total GF G. To 
perform the remaining averaging in the Dyson equation we use the approximation 

G^MG ^ G^MG 

The calculation of M requires further averaging of the product of matrices. We again use 
the prescription of the factorisability there, namely 



M (UmUn) (ImG) (ImG) (ImG) 



However, the quantities UmUn entering into M are averaged here according to 



U^U^ = U2 + (Ui - U2)5^n 

Ui = x^Ul + 2x{l - x)UaUb + (1 - xfUl 

U2 = xUl + (1 - x)Ul (82) 



The averaged value for the self-energy is 



U2 f 

2^ J R{^U^2,i^^)IrnGnrn-a{^l)IrnGrnn-a{^2)ImGrnna{^^) + 



U1-U2 

27r4 



5mn j R{iOi,LU2,LUz)ImGnm-a{^l)IrnGrnn-a{^2)ImGrnna{^z) (83) 



The averaged quantities are periodic, so we can introduce the Fourier transform of them, 
i.e. 

Mmna(a;) = iV"^ ^ M^( A), Cu) e^^{lk{Rm - Rn)) 
k 

and similar formulae for G and G^. Performing the configurational averaging of Dyson 
equation and Fourier transforming the resulting expressions according to the above rules, 
we obtain 

Gr,,{u) = {uj- e{k) - E'^(cu) - M,{k,uj))-^ (84) 

where 

M^{k,Uj) = J R{^l,^2,l^i)N~'^ImGp-q^„{ljJi)ImGq-a{^2) 



pq 



[U2lmGk+pa{^3) + E ImGk+p-a{u;s)] (85) 

9 

The simplest way to obtain the exphcit solution for the self-energyM is to start with 
suitable initial trial solution as it was done for the periodic case (33). For the disordered 
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system, it is reasonable to use as the first iteration approximation the so-called Virtual 
Crystal Approximation(VCA): 

vr 

where for the binary alloy A^Bi^x this approximation read 



£^ = xe^ + (1 - x)el 

Note, that the using of VGA here is by no means the solution of the correlation problem 
in VGA. It is only the using the VGA for the parametrisation of the problem, to start with 
VGA input parameters. After the integration of (83) the final result for the self-energy is 

Mik X _ <E;^q)[^ - njE^n - n{El^,)] + n(E,%^)n(E^-) 

'""^ N'i^ ^ + e;^,-e-^-ei^^ 

{Ui - U2) <E;+q)\^ - n{E~n - n{El^^_^)] + n{El^^_^)n{E--) 
^ uj + E;^^ - E-'^ - El^^_^ 

It must be emphasized that the equations (84) - (85) give the general microscopic self- 
consistent description of inelastic electron-electron scattering in alloy in the spirit of the 
GPA. We take into account the randomness not only through the parameters of the 
Hamiltonian but also in a self-consistent way through the configurational dependence of 
the self-energy operator. 



7 Electron-Lattice Interaction and MTBA 

In order to understand quantitatively the electrical, thermal and superconducting proper- 
ties of metals and their alloys one needs a proper description an electron-lattice interaction 
too A systematic, self-consistent simultaneous treatment of the electron-electron 

and electron-phonon interaction plays an important role in recent studies of strongly corre- 
lated systems. It was argued from the different points of view that in order to understand 
quantitatively the phenomenon of high-temperature superconductivity one needs a proper 
involving of electron-phonon interaction, too |7^- |8^. A lot of theoretical searches for the 



relevant mechanism of high temperature superconductivity deal with the strong electron- 
phonon models. This mechanism is certainly valuable for bismuthate ceramics |Q and for 
fuUerens. Recently a new family of quaternary intermetallic LuNi2B2C compounds 



has been shown to display superconductivity with Tc = 16.6 K for LuNi2B2C, which, 
besides that studies of their physical properties are still in the early stages, suggest that 
electron-phonon coupling is responsible for the superconductivity. 

The natural approach for the description of superconductivity in such type of compounds 
is MTBA [^. The papers |Q,|3^ contain a self-consistent microscopic theory 



of the normal and superconducting properties of transition metals and strongly disor- 
dered binary alloys in the framework of Hubbard Model (11) and random Hubbard model 
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(73). It is worthy to emphasize that in paper |34| a very detailed microscopic theory of 
strong couphng superconductivity in highly disordered transition metals alloys has been 
developed on the basis of IGF method within MTBA reformulated approach [Q. The 
Eliashberg-type strong couphng equations for highly disordered alloys has been derived. 
It was shown that the electron-phonon Spectral Function in alloy is modified strongly. An 

H 



interesting discussion 

(c.f. PI- m). 
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clarified many uncertainties in this important issue 



8 Other Applications of the IGFs Method 



Another important application of IGF method is related with the investigation of non-local 

A comparative 



correlations and quasiparticle interactions in Anderson model |2^ 
study of real-many body dynamics of single-impurity, two-impurity and periodic Anderson 
model, especially for strong but finite Coulomb correlation, when perturbation expansion 
in U does not work (c.f. |^) has permitted to characterize the true quasiparticle excita- 
tions and the role of magnetic correlations. It was shown that the physics of two-impurity 
Anderson model can be understood in terms of competition between of itinerant motion of 
carriers and magnetic correlations of the RKKY nature. The correct functional two-pole 
structure of the propagator has been found for the strongly correlated case. This issue is 
still very controversial 
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and the additional efforts must be applied in this field. 
The application of the IGF method to the theory of magnetic semiconductors was very 
succsessfull |28]. As a remarkable results of our approach let me mention the fi- 



nite temperature generalization of the Shastry-Mattis theory for magnetic polaron |28 



which clarified greatly the true nature of the carrier in magnetic semiconductors. There 



are some analogy of the Kondo-lattice type of model in ||2^ with the Kondo-Heisenberg 
model of copper oxides, however the physics are different. There is a dense system of 
spins interacting with smaller concentration of holes in HTSC. The application of IGF 
method to spin-fermion model | p.OU|| has been done by using the theory of Heisenberg 
antiferromagnet and allows one to consider carefully the true nature of the carriers 
in Cu02 planes. 

And finally, the new interesting application of the IGF method for consideration of dy- 
namics of quasiparticles and dynamical conductivity of single electron resonant tunneling 



systems has been done recently in papers ||101|| , ||102|| (c.f. [|103|| ). This reformulation of 



IGF method has much in common with the approach of paper 104 



9 Conclusions 

In the present paper we have formulated the theory of the correlation effects using the 
ideas of the quantum field theory for the interacting electron system on a lattice. The main 
achievement of this formulation is the derivation of the Dyson equation for two-time ther- 
modynamic retarded Green's Functions instead of causal ones. Such a formulation permit 
to use the convenient analytical properties of retarded and advanced GF and advantage 
of using the formal solution of the Dyson equation, which, in spite required approxima- 
tions for the self-energy, provide the correct functional structure of the single-electron GF. 
This strong point of our approach do not give the possibility of direct application of it to 
the calculation of the two-particle GFs. In this paper we have considered in details the 
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idealized single-band Hubbard model, which is one of the simplest (in the sense of formu- 
lation, but not solution) and most popular model of correlated lattice fermions. We have 
presented here the novel method of calculation of the quasiparticle spectra for this model, 
as the most representative example. We hope that this explanation have been done with 
sufficient details to bring out their scope and power, since we believe that such techniques 
will have application to a variety of many-body systems with complicated spectrum and 
strong interaction, as it was shown in Section 7. 

In summary, with using IGF method we were able to obtain the closed self- consistent 
set of equations determining the electron GF and self-energy. These equations define the 
renormalization coefficient of the one-electron GF0, defined for a point (k, uj = e{k)): 

The renormalization coefficient (87) is an one of the most important notion for the char- 
acterization of the single-particle behaviour of the quasiparticle excitations in correlated 
many-body systems. For the Hubbard model, these equations give the general microscopic 
description of correlation effects for both the weak and strong Coulomb correlation, de- 
termining of the complete interpolation solution of the Hubbard model. Moreover, this 
approach gives the workable scheme for the definition of the relevant Generalized Mean 
Fields written in terms of appropriate correlators. The most important conclusion to 
be drawn from of the present consideration is that the GMF for the case of strong 
Coulomb interaction have a quite non-trivial structure and cannot be reduced to the 
mean-density functional. This last statement resemble very much the situation with the 
strongly nonequilibrium system, where the single-particle distribution function only not 
enough to describe the essence of the strongly nonequilibrium state and more compli- 
cated correlation functions must be taken into account, in accordance with general ideas 
of Bogolubov and Mori-Zwanzig. The IGF method is intimately related to the projection 
method in this sense, which express the idea of a "reduced description" of the system in 
the most general form. This line of consideration are very promising for developing the 
complete and self-contained theory of the strongly interacting many-body systems. 
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